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Abstract 

We present a calculation of the cross section for the process e + e~ — > W^H^ 
in the Minimal Supersymmetric Standard Model (MSSM) and the Two Higgs Dou- 
blet Model (THDM). We study the basic features of the MSSM prediction for some 
distinctive parameter scenarios. We find large effects from virtual squarks for sce- 
narios with large mixing in the stop sector which can lead to a cross section vastly 
different from a THDM with identical Higgs-sector parameters. We investigate this 
interesting behaviour in more detail by thoroughly scanning the MSSM parameter 
space for regions of large cross section. For a charged Higgs boson too heavy to be 
pair-produced at such a machine, it turns out that a large MSSM cross section with 
a good chance of observation is linked to a squark mass scale below 600 GeV and a 
considerable amount of mixing in either the stop and sbottom sector. 



1 Introduction 



The discovery potential of the CERN Large Hadron Collider (LHC) should be sufficient 
to resolve the issue of the existence of a neutral Higgs boson pp. An electron-positron 
collider would, however, serve as an ideal tool to measure the properties of neutral Higgs 
bosons very precisely [2]. To discover a charged Higgs boson, H ± , on the other hand, is 
much harder at the LHC [TJ [3] , especially if the H is substantially heavier than the top 
quark. The reason for this is the dominant decay of the H bosons into heavy quarks (tb or 
ib) which leads to hadronic signatures that are hard to distinguish from QCD background 
events and which forces one either to rely on the less probable decay H — ► tu or to 
accumulate a lot of statistics during several years of running at high luminosity in order 
to claim discovery Actually, the strong dependence of the production cross sections on 
the ratio of vacuum expectation values in the Higgs sector, 1)2 /v\ = tan f3, makes it almost 
impossible to discover the charged Higgs boson at the LHC in a wedge-shaped region 
of intermediate tan f3- values in the tan/3-m#± plane if it is heavier than the top-quark 
[HE]. This non-discovery range includes 4 < tan (3 < 12 at m H ± = 250 GeV and widens 
continuously to e.g. 2 < tan/? < 40 at m H ± = 650 GeV [3]. 

At an e + e~ collider with a centre-of-mass energy y/s, the main production process 
for charged Higgs bosons is pair production (e + e~ — > H + H~) which is mediated mainly 
via photon- and Z-exchange in the s-channel. The pair-production cross section is almost 
independent of tan f3 and, owing to the colourless initial state, Higgs-boson signatures from 
decays into heavy quarks are not obscured by enormous amounts of QCD background 
events, like at the LHC. Therefore, detection of charged Higgs bosons via pair production 
should be possible in the whole m#±-tan (3 plane up to a mass limit somewhat below yfs/2. 
If the collider energy is not sufficient for pair production (i.e. \fs < 2m#±), the only way 
to extend the search for the charged Higgs boson to higher masses is to try to observe 
its single production in association with lighter particles. The most relevant processes 
which have been studied in the literature are the tree-level processes e + e~ — > H tu [U [5] 
and e + e~ — > H ± tb P, [5], and the loop-induced processes e + e~ — ► H ± eu El [5] and 
e+e~ -> W ± H T P EH dH [El Q31 [5]. Most of the studies are done in the framework of 
a Two Higgs Doublet Model (THDM) with parameters fixed at MSSM values, except for 
Refs. pH H2] which study W ± H T production in the general THDM, Refs. 0, [13] which 
study the same process in the full MSSM, and Ref. which includes virtual sfermion 
contributions for the cross section prediction of the H ± eu e final state. 

If the collider energy does not allow for pair production and yfs > raw + rn H ± , the 
W ± H T production process can become the dominant source of charged Higgs bosons, 
depending on tan (3. In this the process is loop-induced, the signal rate is much 

smaller than the typical pair-production rate. Even if H + H~ production is the main 
production process at an e + e~ collider, it would still be very rewarding to measure the 
W ± H T production cross section as well. While the H + H~ cross section is determined at 
tree-level by the gauge couplings to photon and Z-boson, the W ± H T cross section is loop- 
induced and thus depends already at leading order on the virtual particles in the loops. 
Therefore, with a cross section measurement of the latter process one would have access 
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to information about the underlying model. 

This paper presents a calculation of the MSSM prediction for the process e + e~ — > 
W + H~. By the time this work was finalized, another calculation of the same prediction 
appeared [9] . We carried out a detailed comparison [TJ] of our results with those of Ref . [9] 
which lead to agreement^. In the framework of the Two Higgs Doublet Model (THDM), 
the cross section prediction for e + e~ — > W + H~ is well-known [TO] [TT] [T2] and we could 
reproduce the results of Ref. [12] in particular. 

In Section [2] we present the process e + e~ — > W + H~ with the contributing Feynman 
diagrams in the MSSM and some details about our calculation. Section [3] presents our 
numerical results in two stages. In stage one, we examine the basic features of the W ± H T 
cross section for a 500 GeV and 1 TeV e + e~ collider for two distinct MSSM parameter 
scenarios and compare the results also to predictions of a THDM with the same Higgs sector 
parameters. This enables us to exemplify the range of influence of virtual superpartners on 
the cross section. The effect of polarized electron and positron beams on the cross section 
is discussed briefly followed by a demonstration of the decoupling of superpartners with 
increasing supersymmetry-breaking scale. In stage two, we present results of a thorough 
scan of the MSSM parameter space for regions of large, possibly observable, cross section 
and try to understand for which parameter scenarios they arise. Section H] contains our 
conclusions. Some results of this calculation were already reported in [T3] . 



2 e+e" W+H- in the MSSM 



2.1 Kinematics 

We study the reaction 

e + (k, a) + e~(k, a) — > W + (p, A) + H~ (p) , 

where k and k denote the momenta of the initial-state electron and positron, p and p 
the momenta of the final-state gauge boson W + and Higgs boson H~ . Additionally, the 
electron, positron, and IV-boson are characterized by their spin polarization a, er(= ±|) 
and A(= 0, ±1). Neglecting the electron mass, the kinematical invariants s = (k + k) 2 , 
t = (k — p) 2 and u — (k — p) 2 fulfil the relation 

s + t + u = mly + m 2 H ± . (1) 

We assume unpolarized electron and positron beams. Thus, the differential cross section 
summed over spin polarizations of the final-state H^-boson reads 



da 1 



dt 16vrs 2 ^ 4 

A=0,±1 a,a=±\ 



1 As a result, the formulae in [9] are now checked by an independent calculation. According to [14] the 
results agree, if in Eq. (C14) of Ref. [H] the tensor coefficient D 2 3 in the coefficient of A^g^ is replaced 
by -ID,;. 
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with the helicity amplitudes M. a a\ of the process. The total cross section is evaluated by 
numerical integration over the kinematically allowed t-range: 

/■tmax(s) ^ 

a e+e -^ w+H -{s) = / dt—(s,t). (3) 

2.2 Feynman graphs 

In the MSSM and THDM all leptons couple gauge-invariantly to one of the two Higgs 
doublets and therefore also to the physical Higgs bosons. Thus, there is a non-vanishing 
tree-level amplitude for the process under study. The Feynman graphs on tree-level consist 
of three graphs with s-channel exchange of neutral Higgs bosons (h°, H°, A ) and one with 
a t-channel exchange of a neutrino. The Yukawa couplings which appear in all tree-level 
graphs are oc m e /my/ pa 6 ■ 10~ 6 and suppress the tree-level contribution strongly. The 
tree-level amplitude can thus safely be neglected and consequently there are no photon 
bremsstrahlung corrections at leading order. 

In our calculation we take into account all one-loop contributions to the amplitude 
which do not vanish in the limit m e = 0. For this reason, Feynman graphs with insertion 
of a s-channel Z—A mixing self-energy, or a t-channel neutrino self-energy, or radiative 
corrections to e + e~{h°, H°, A } Yukawa couplings or to the e ± v e H T Yukawa coupling need 
not be considered. There remain Feynman graphs with insertions of a W~-H~ or G~-H~ 
mixing self-energy at the external leg of the outgoing charged Higgs boson (see Fig. [I]), 
graphs containing the loop-induced r )W + H~ and ZW + H~ vertices^ and box-type graphs 
(see Fig. [2J and EJ). The amplitude divides into Feynman graphs which contain either solely 
THDM particles (Fig. [Ib, c upper lines and Fig. [2}, or solely superpartners (Fig. [Tb,c lower 
lines and Fig. [3D in the loop. 

2.3 Calculation 

Although the tree-level contribution vanishes in the the limit of vanishing electron mass, 
which we consider, the need for renormalization arises at one-loop. There are divergent 
vector-scalar and scalar-scalar mixing propagators (Fig. [T]) and vertex-type graphs (Fig. [2] 
and [3}[1 We use the on-shell renormalization scheme of Ref . [16] , the application of which 
to W ± -H T self-energies is discussed in detail in Ref. [17]. In the MSSM no new types 
of divergent vertex functions occur. Merely, additional loop contributions with virtual 
superpartners add to the existing types of vertex functions. Therefore, we can take over 
the renormalization conditions from Ref. [UJ directly, which are briefly summarized in the 
following. 

In order to generate the counter-terms, which are needed to renormalize the amplitude 
under study, it is sufficient to introduce field renormalization constants for the MSSM 

2 Note that there is no tree-level ZW ± H^ coupling in the MSSM or THDM [T5] . 

3 In principle, also divergent tadpole insertions appear which we do not display, because we will choose 
a renormalization condition such that those contributions vanish (see Eqs. © and (jTUJ) ) . 
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(c) G + -H self-energy 

Figure 1: Feynman graphs for the process e + e~ — > W + H~ with self-energy insertions. 
Each combination of particle labels corresponds to a separate Feynman graph. 
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Figure 2: The THDM subset of Feynman graphs of vertex and box type for the process 
e + e~ — > W + H~ . Each combination of particle labels corresponds to a separate Feynman 
graph. 
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;ure 3: The superpartner-loop subset of of vertex and box type for the process e + e 
V H~ . Each combination of particle labels corresponds to a separate Feynman graph. 
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Higgs doublets Hi, H 2 and renormalization constants for their vacuum expectation values 
Vi, v 2 : 

Hi -> y/Z^Hi , Vi -> y/Z^(Vi - <J«i) , z = 1, 2 . (4) 



Expanding the renormalization constants in the MSSM Lagrangian to one-loop order, 
Z H . = 1 + SZff' and 8vi = 5v[ 
corresponding Feynman rules: 



Zjj. = 1 + 5Zft and 8vi = 5v[ , generates the counter-term interactions needed and the 



rcrlH^W*^)] = i m 2 w 6Z HW , (5) 

TcTijuW^H^] = -iem w g^ 5Z HW , (6) 



T UT [Z ll W^H^}=iem w —g llv 5Z HW (7) 

Cm 



with 

K7 ■ n nf^i 5v 2 6Z H2 -6Z Hl 
oZhw — sm p cos p h 



v 1 v 2 2 

where Z, H^ und 7 denote the electroweak gauge bosons and the photon, and the 
momentum of the W ± boson, chosen as incoming. In the on-shell scheme the following 
renormalization conditions are posed: 

• Renormalized tadpole graphs vanish, i.e. 

th° = th° + 8th° = , (9) 
t H o = t H o + 5t H o = . (10) 

This guarantees that the parameters Vi, v 2 in the renormalized Lagrangian describe 
the minimum of the Higgs potential at one-loop order. 

• Real charged Higgs bosons H ± do not mix with longitudinally polarized W ± bosons, 
i.e. the real part of the renormalized H ± -W T mixing self-energjfl, 

T>Hw(k 2 ) = Y> HW (k 2 ) - mlySZnw , (11) 

vanishes if the momentum k of H ± is on mass-shell: 

Re± HW (k 2 )\ k2=m2 =0. (12) 

The condition ffl2l) fixes the renormalization constant SZhw'- 

5Z HW = -^-ReY; HW (m 2 H ±) . (13) 



in 



w 



The renormalized H ± -W T mixing self-energy is denned as the coefficient of — i;^— of the amputated 
renormalized H^-W^ propagator. 
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Figure 4: Counter-term diagrams for e + e — ► W + H . 

Because of the BRS symmetry of the renormalized Lagrangian the renormalization of the 
divergent H^-G^ mixing self-energy is connected with the H ± -W T mixing self-energy 
through a Slavnov-Taylor identity [IT] : 

k 2 t HW (k 2 ) - m 2 w ± HG (k 2 ) = . (14) 

As a consequence, the real part of the renormalized H ± -G T mixing self-energy, 

t HG {k 2 ) = T, HG {k 2 ) - k 2 5Z HG , (15) 

also vanishes for k 2 = m 2 H± : 

Re£ HG (k 2 )\ k2=ml± =0. (16) 

The Feynman rule for the corresponding counter-term interaction thus reads: 

T C T[H^G ± (k)]=ik 2 6Z HG (17) 

with 5Z HG = -ReJ: HG (m 2 H± )/m 2 H± . 

The calculation of the amplitude has been performed in 't Hooft-Feynman gauge using 
Constrained Differential Renormalization [18J with the help of the computer programs 
FeynArts and FormCalc [19J. To that end, the counter-term definitions and Feynman rules 
of counter-term interactions which were needed were added to the MSSM model file for 
FeynArts. 

3 Results 

Our presentation of results is meant to be exemplary. Given the possible time-line for the 
realization of a high-energy e + e~ collider like the International Linear Collier (LLC) and 
given that it is not yet known what the LHC will reveal about this new energy frontier, 
there is no point in attempting a comprehensive presentation. The predictions for W ± H T 
production in the framework of the general THDM are well known [TO f fl~Tj [I~2] . Our concern 
is to demonstrate possible distinctive differences between the MSSM and the THDM. To 
that end, we compare the prediction of our complete one-loop calculation in the MSSM 
with the THDM prediction with the Higgs-sector parameters restricted to MSSM values 
(sTHDM). The latter model corresponds to the MSSM with decoupled superpartners. 
Thus, one achieves insight in the contribution of superpartner loop graphs to the cross 
section. 
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3.1 MSSM parameter scenarios 



In order to demonstrate possible differences between the sTHDM and the full MSSM, we 
are interested in scenarios which potentially show large effects from virtual superpartners. 
Hence, the MSSM parameter scenarios for which we present numerical results have a rather 
low squark-mass scale. In the MSSM, large mixing in the stop sector is known to induce 
large radiative corrections in the Higgs sector. Therefore, we investigate two types of 
scenarios, one with large mixing in the stop sector and the other one with no mixing. 

In order to specify the relevant MSSM parameters, let us recall the stop mass matrix. 
For real values of the supersymmetric Higgsino mass parameter \x and the soft-breaking 
parameter A t the stop mass matrix reads in the L R basis 



where X t = A t — /icot/3, Q t = 2/3, and s w = sin9 w . Mq and Mq are the soft-breaking 
masses for squark isospin-doublets and -singlets respectively. The two MSSM parameter 
scenarios are specified as follows: 

Large- stop-mixing scenario: We set all sfermion mass parameters (Mq, Mjy, Mg, Mi, Mg) 
equal to one common sfermion mass scale = 350 GeV, the stop mixing parameter 
X t = —806 GeV, fi = 300 GeV, the gaugino mass parameters M 2 = 1000 GeV and 
M 1 = (5s 2 w /3c 2 w )M 2 = 476.26 GeV (GUT relation). 

The main features of this scenario are: 

• The mass of the light stop is of order 100 GeV and there is a large mass splitting 
between the two stop mass eigenstates. 

• The stop mixing angle is approximately 45°. Thus, maximal mixing occurs, i.e. 
all entries in the stop mixing matrix have approximately the same magnitude. 

• The mass of the lightest MSSM Higgs boson mh is almost maximal with respect 
to X t (because \X t \ w 2M f - [20]), e.g. for m H ± = 350 GeV and tan/3 = 30 we 
have m h = 120.06 GeV. 

• This scenario allows for the production of stop pairs at the energy scales con- 
sidered in the following sections and the decay of the charged Higgs into stop 
and sbottom will be kinematically allowed for m#± > 450 GeV. 

No- stop-mixing scenario: Here, we choose the same values of the MSSM parameters as 
above, except for X t which is set to zero. 

The main features of this scenario are: 

• Both stops have masses around 400 GeV and all other sfermions have masses 



M~ 




+ m\ + m|(| - Q t s 2 w ) cos 2/3 



m t X t 



m t Xt 

M} + m\ + m 2 z Q t s 2 w cos 2/3 




around 350 GeV. 
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• The stop mixing angle is zero. Thus, no mixing occurs, i.e. the L R basis for 
stops coincides with the mass basis. 

• The mass of the lightest MSSM Higgs boson is rather low, e.g. = 103.02 GeV 

for m H ± = 350 GeV and tan/5 = 30. 

For the values of and tan/3 chosen for our discussion, both scenarios respect all 

superpartner mass bounds from direct search results and fulfil indirect constraints on the 
parameters coming from requiring vacuum stability and from experimental bounds on the 
electroweak precision observables. Also the lower bound on the lightest neutral MSSM 
Higgs boson is respected by both scenarios, except for the case tan (3 = 3 in the no-mixing 
scenario. Yet, for the sake of comparison, we also show this case among our numerical 
results. 

3.2 Cross-section for the parameter scenarios 

There are two major motivations for studying the MSSM prediction for e + e~ — > W ± H T 
cross sections at a future ILC Firstly, it is important to know the expected event rates at 
the ILC either to confirm the discovery of charged Higgs bosons at the LHC or to assess 
the H discovery potential of the ILC beyond the pair-production limit of \fs/2. Secondly, 
the observation of this process provides the opportunity to glean some information about 
the underlying model because of its potential sensitivity to virtual superpartners. 

We present results for two collider energies, y/s = 500 GeV and 1000 GeV. For each 
energy we show results for the the large-stop-mixing and no-stop-mixing scenario for dif- 
ferent tan /3- values in and around the "wedge region" and for charged Higgs bosons masses 
below, above and right at the pair production limit of 250 and 500 GeV, respectively. 

3.2.1 Large-stop-mixing scenario 

In Figs. [5^, and [6^ the e + e~ —>■ W ± H T cross sections for a collider energy of 500 GeV 
are displayed for the full MSSM and the corresponding sTHDM. Fig. [5^ shows the m^±- 
dependence for tan/3 = 3, 15, and 30. The spike in all predictions for m H ± = y/p* m 
180 GeV m mj + is due to a threshold effect in the vertex graphs with virtual top- 
and bottom-quarks (see Fig. [2]). Similar spikes appear in the full MSSM prediction in 
places where m H ± equals + m~ b .(i,j = 1,2) and are due to threshold effects in the 
Feynman graphs with virtual stops and sbottoms, e.g. for tan/3 = 30 where « 
390 GeV ~ mi x + m~ b . Fig. shows the tan /3-dependence for m H ± = 190 GeV, 250 GeV, 
and 350 GeV. The MSSM cross section is generally about two orders of magnitude larger 
than the sTHDM in the large-stop-mixing case for large values of tan /?. 

For */s = 1000 GeV Figs. [5b and [6b reveal the same general behaviour. Yet, due to the 
higher collider energy, more stop-sbottom thresholds are within the accessible mass range 
for the charged Higgs boson. 

The enhancement with respect to the sTHDM case is due to the threshold effects 
combined with enhanced couplings of the charged Higgs boson to third-generation squarks. 
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In the large-stop-mixing scenario, the Feynman graphs with virtual squarks of the third 
generation comprise the dominant contribution to the cross section if tan (5 is not too small. 
In this scenario each entry of the rotation matrix which transforms the stop fields into the 
mass eigenbasis is of similar magnitude. Thus, all potentially large terms in the couplings 
of third-generation squarks to Higgs bosons^l are weighted by a factor of order 1, especially 
the term proportional to m t \i. 

3.2.2 No-stop-mixing scenario 

Figs. [5b and [6b show the predicted cross section for a 500 GeV collider in this scenario. 
Superpartner effects are small in this case and the MSSM as well as the sTHDM predict 
almost the same cross section. Without mixing in the stop sector the two factors which 
caused the large enhancement in the large-mixing case are absent. Firstly, there are no 
enhancement effects by thresholds in stop-sbottom loops. Secondly, the mixing angle in 
the stop sector is exactly zero which leads to the absence of terms proportional to m t fi in 
the t\biH -couplings. 

For a/s = 1000 GeV the situation changes (see Figs. [5bl and EH). The maximally 
accessible charged-Higgs-boson mass is higher than all thresholds in stop-sbottom loop 
graphs (730 GeV < m^. +m^. < 760 GeV). Again, there are threshold enhancements of the 
MSSM prediction with respect to the sTHDM of up to one order of magnitude for large 
tan f3. 

3.2.3 Effect of polarized electrons and positrons 

The cross section of e + e~ — ► W ± H T depends strongly on the polarization of the incoming 
electrons and positrons. Specifically, the cross section varies roughly by one order of 
magnitude between the optimal situation where all electrons have helicity —1 and all 
positrons have helicity +1, i.e. P(e~) = —1 and P(e + ) = +1, and the opposite situation 
where P(e~~) = +1 and P(e + ) = —1. 

In order to illustrate the effect of polarized beams at a future LC, we include in Figs. [5b. 
to [5bi the cross section prediction for the optimal polarization for tan/3 = 15. Likewise, 
we include in Figs. [6b to [6bl the cross section prediction for the optimal polarization for 
m H ± — "v/s/2. Optimal polarization leads to an increase in cross section varying between 
a factor of 2 up to 3.8 depending on m H ± and tan/? for both collider energies. 

3.2.4 Decoupling of superpartners 

We set all soft-breaking parameters which appear in our^expression for the MSSM cross 
section equal to a single supersymmetry breaking scale M, 

M = Mf = A b = A t = Mi = M 2 = n , (19) 

5 A full list of the relevant couplings in our notation can be found in [21] 
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and study the behaviour of the MSSM cross section with rising M compared to the sTHDM 
prediction. This will allow us to demonstrate the decoupling of the virtual superpartner 
contributions to the cross section. 

Fig. [7h shows results for a 500 GeV collider and m#± = 190, 250, and 350 GeV with 
M varied between 350 and 2000 GeV. The MSSM results drop significantly once all super- 
partners get heavy enough for threshold effects to disappear. A similar behaviour is shown 
in Fig. [7b for a 1000 GeV collider and m H ± = 350, 500, and 650 GeV where M is varied up 
to 4TeV. Notably, even at a soft-breaking scale of ITeV, the MSSM surpasses the sTHDM 
prediction by a factor of 1.4, 1.7, and 2.2 for these three mass values, respectively. 

3.3 Parameter scan 

We demonstrated the possibility of strong enhancement of the MSSM cross section predic- 
tion by virtual superpartner effects for the large-stop-mixing scenario. This might predict 
W ± H T production to be observable in such MSSM scenarios unlike in the corresponding 
sTHDM scenario. We investigate this possible behaviour more systematically by scanning 
over the relevant MSSM parameter space. 

The detectors at the LHC will be capable of finding an MSSM charged Higgs boson 
only in certain regions in the m^o-tan (3 plane. Likewise, the ILC can detect charged Higgs 
bosons easily only if pair production is kinematically possible. Thus, as far as W ± H T 
production at the ILC is concerned, there are two points of departure: either we will know 
to some extent what the mass of the charged Higgs boson is or we will not. In the first 
case it makes sense to combine this knowledge with the knowledge of the cross section 
prediction in the MSSM and adjust the collider energy to the maximum of the expected 
cross section as a function of \/s. In the second case it is difficult to say what the best 
strategy of choosing \fs will be for searching for the charged Higgs boson. A reasonable 
assumption is however that a significant fraction of the data during ILC operation will be 
collected at the highest collider energy. 

With these two possibilities in mind, we devised two parameter scans. The first scan 
assumes that the charged Higgs mass is known and that it makes sense to adjust the collider 
energy to maximize the cross section. Furthermore, it assumes that one has also acquired 
some knowledge about the value of tan f3. We thus fix m#± and tan /3 in our scan and 
include y/s in the list of scanned parameters. The second scan fixes y/s at the maximal 
collider energy but includes and tan f3 in the list of scanned parameters. For both 

scans we assume an ILC with a maximal collision energy of 1 TeV. 

In general, we are interested in regions of parameter space where the cross section is 
large such that it might still be detectable at the ILC. In order to increase the number 
of scanned parameter points in the region of large cross section, we use the "adaptive 
scanning" algorithm described in more detail in Ref. [22]. Essentially, that means that 
we are using the importance sampling algorithm VEGAS [23] to evaluate the integral of 
the cross section over the relevant MSSM parameters. However, we are not interested in 
the value of the meaningless integral but store the sampled points which automatically 
accumulate in the region of large cross section. This method allows to perform scans of 
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higher-dimensional parameter spaces with emphasis on certain features of the results much 
more efficiently than a grid of points or purely random sampling can do. More precisely, 
the quantity we are integrating over in our scans is identical to the cross section only at 
'allowed' points according to the exclusion limits defined below, and zero otherwise. The 
enrichment of points thus focuses also on the allowed region. We use the following exclusion 
limits: 

m- u > 95.7 GeV , m~ b . > 89 GeV , m~^~ > 150 GeV , 

m^o ^ 46 GeV , m^± ^ 94 GeV , M i > 73 GeV , (20) 

\Ap\ < 0.0025, 

where Ap is the dominant squark contribution to the electroweak rho-parameter. The 
mass bounds on superpartners and on Ap are according to Ref. [24J. 

Apart from these exclusion limits, we take into account two more major constraints on 
the MSSM parameter space. 

Firstly, we calculate for each parameter point the MSSM predictions for m^o and 
<r(e + e~ — > h°Z) x BR(/i° — > bb) and exclude it if the m^o-dependent LEP-bound on 
a x BR is violated (according to Table 14(b) of Ref. [25]). We use FeynHiggs 2.5.1 [26] for 
calculating the m^o prediction and allow for a theoretical uncertainty of 3 GeV. 

Secondly, we calculate the leading-order MSSM prediction for the branching ratio 
BR(B — >• X s 7) [27J and exclude parameter points if the prediction falls outside of the 
range (3.55 ± 1.71) • 10~ 4 . This range is determined by using the experimental central 
value [2B] and adding up the experimental 3a interval (~ 10~ 4 ) and an estimate of the 
independent theoretical uncertainty (0.71 ■ 10~ 4 ). The latter estimate (20%) is guided by 
the discussion of theoretical uncertainties for the SM prediction |29j. 

3.3.1 Scan 1: m#± and tan/3 known 

The MSSM input parameters are scanned over the following region: 

Mj, Mi, M 2 = 10 . . . 2000 GeV , 
/i, A t , A b = -4000 . . . 4000 GeV , 

y/s = 500... 1000 GeV, (21) 
m H± = 500 GeV , 
tan/3 = 3, 15, 30. 

The tan f3- values respectively lie roughly at the lower, middle, and upper end of the wedge- 
region. Examining the results by studying projections of the cross section on all one- and 
two-dimensional subspaces of the scanned parameter space we find that the gaugino mass 
parameters Mi and M 2 have negligible influence on the cross section, at least in all regions 
where the cross section is above 10~ 3 fb. Because of their importance we also study the 
dependence of the cross section on the squark mixing parameters X t , X\, normalized to the 
sfermion mass scale Mj, i.e. X t = (A t — pj tan /3)/Mj and X b = (A b — /itan/3)/M?. In 
Fig. [S]we show the most interesting two-dimensional projections [3Ti] . 
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For tan (3 = 3 the most noticeable feature, shown in Figs. [Hb, [Hb and [Hb, is that the 
highest values for the H production cross section lie typically between 10~ 3 fb and 
1CT 2 fb and are reached everywhere in the allowed parameter region. In our scan we find 
a few cases with higher cross section, clearly visible as the scattered yellow dots on a cyan 
background. Yet, we find no scenario with a cross section above 0.1 fb. The well-known 
fact that for small tan f3 the constraints, especially the LEP Higgs-mass bound, require a 
fair amount of stop mixing (i.e. \X t \ > 1) to be present, is clearly seen in all three figures. 

Also owing to the LEP bound, the lowest common squark mass scale Mj for which we 

find allowed scenarios is about 400 GeV for positive X t and 600 GeV for negative X t . It 
is for rather low Mj, below about 650 GeV, that the few scenarios with a cross section 
above 10" 2 fb appear. It turns out that having Mj < 650 GeV is a rather generic necessary 
requirement for scenarios with cross sections above 10~ 2 fb, as will become more clear in 
the following. 

The projection of the parameter points on the X t —\i plane in Fig. [Hb reveals that almost 
all of the scenarios we find in our scan show a strict correlation between the sign of X t and 
the sign of fi, the former being positive if the latter is negative and vice versa. 

For tan (3 = 15 (Figs. Ebl, Eb andEF) the results of the scan change drastically compared 
to the tan/3 = 3 case. Firstly, the allowed region in the Mf-X t plane (see Fig. [HJl) is 
much larger, now also including scenarios without stop mixing and scenarios with Mj 
below 300 GeV. Secondly, the highest cross section values are considerably larger, of the 
order of 10 fb. Interestingly, the projection of the parameter points on the Mj-X t plane 
(Fig. [Hbl) shows that the scenarios with a (large) cross section above 1 fb are confined to 
a very specific region with between 250 GeV and 600 GeV, and X t between —2.5 and 
— 1 simultaneously. 

The projection on the X t -Xb plane (Fig. [Hb) shows that all large cross section scenarios 
also have a significantly nonzero and negative value of X^. Quantitatively, it turns out 
that those scenarios typically have the two mass scales mtX t and rribXb of the same order 
of magnitude. Furthermore, from the projection on the X t -[i plane (Fig. [HJ) we see that 
those scenarios also have a positive and large value of /i, typically between 2 TeV and 4 TeV 
(the limit of our scan range). Looking into the Xb-Ab plane (not depicted) we find large 
cross section scenarios for almost any value of Ab. 

The scan results for tan /3 = 30 (Figs. [Hg, [HJi and [HJ) look very similar to the previous 
case. Again, in order to obtain a cross section above 0.1 fb it appears that one needs 
a low sfermion mass scale Mi between roughly 250 GeV and 600 GeV and a significant 

amount of mixing in the stop and sbottom sector with both X t and Xb negative. One 
important change compared to the previous case is that all large cross section scenarios 
have a rather large value of | A b \ > 1 TeV (see Fig. [HJ) while the values of fi are constrained 
to 0.8 TeV < /i < 2.2 TeV (result not depicted). 
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3.3.2 Scan 2: m#± and tan/3 unknown 



The MSSM input parameters are scanned over the following region: 



Mi = M 2 



M f 
H,A t ,A b 



m H ± 
tan (3 



10...2000GeV, 
-4000 . . . 4000 GeV 
500...920GeV, 
1 . . . 50 , 
500 GeV , 
ITeV. 



(22) 



The less influential gaugino-mass parameters Mi, M 2 have been fixed in order to obtain a 
more thorough scan over the relevant MSSM parameters. 

Fig. [9^l shows the obtained cross section values projected on the m#±-axis. Overall, the 
values drop with rising m#± due to decreasing phase space. There are no scenarios with 
a > 0.1 fb for m#± > 750 GeV and no scenarios with a > 0.01 fb for m#± > 875 GeV. The 
projection on the tan /3-axis in Fig. [Hb shows that there exist allowed parameter scenarios 
with cross sections even above 0.1 fb over the whole scanned range. Fig. [9b shows the very 
distinct dependence of the cross section on the common sfermion mass scale Mj which 
survives even after the inclusion of m H ± and tan/3 in the scan. Evidently, cross sections 
significantly above 0.1 fb are only possible if 275 GeV < Mj < 700 GeV. A detection of 
squarks with masses in this range at the LHC would raise hopes that W ± H T production 
would be observable at the ILC with an integrated luminosity of the order of 1000 events/fb. 

As in Scan 1 it turns out that the mixing parameters in the third generation squark 
sector are important for the understanding of the large cross section regions. Fig. [TUk shows 
that the largest cross section values lie roughly on two broad stripes with 1 < \X t \ < 2.5 in 
the X t -Xb plane leaving the central part of the plane with values below 0.01 fb. We find 
cross section values above 0.1 fb only for scenarios which also have a considerable amount 
of sbottom mixing, typically Xt < —50. 

At this point, it is instructive to compare our scan with a variant version where the LEP 
Higgs constraint on the MSSM parameters is implemented in its most simple way, which 
is often done in the literature. This simple implementation only makes use of the overall, 
scenario independent, bound on m h o (e.g. 89.8 GeV [21]), discarding all scenarios which 
predict m^o to lie below it. Obviously, this prescription is overly conservative, keeping e.g. 
all SM-like scenarios with masses above 89.8 GeV although in this case a bound close to 
the SM limit of about 114.4 GeV [31] applies. 

Fig. [TUb shows the results for Scan 2 using the simple implementation of the LEP Higgs 
constraint. The result looks quite different to the detailed implementation of the constraint 
in Fig. [TUk . The largest cross section values lie roughly scattered around a half-ring in the 
X t -Xb plane. In this less restricted version of Scan 2 a certain amount of squark mixing, 
either in the stop or sbottom sector is needed in order to obtain cross sections above 0.1 fb. 
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Motivated by this observation we study the dimensionless quantity 



X t b '■ 




2 



m h X h 



) 



2 



(23) 



as a discriminating variable for our original Scan 2. This quantity is just the off-diagonal 
entries in the stop and sbottom mass matrices, m t X t and rribX b respectively, added in 
quadrature and normalized to the common sfermion mass scale Ms. Fig. [TUb shows the 

projection in the Mj-Xtb plane. Evidently, cross sections above 0.1 fb require X t b > 0.8 

and centre roughly around Xtb ~ 1- 

4 Conclusions 

The production of a charged Higgs boson H ± in association with a an electroweak boson 
W T in e + e~ collisions is a loop-induced process. At the ILC, this process is particularly 
important for the charged- Higgs-boson detection if pair production is kinematically for- 
bidden. Furthermore, one would gain valuable information about the underlying model 
by observing this process at the ILC. We calculated the cross section for this process at 
one-loop order in the framework of the MSSM and THDM and investigated the predictions 
of both models at an ILC with 500 GeV and 1000 GeV centre-of-mass energy. The MSSM 
scenarios with large stop mixing and low sfermion mass scale, for which we showed an exam- 
ple, can give rise to a cross section which differs from a THDM with identical Higgs sector 
parameters by two orders of magnitude. Using polarized electron and positron beams can 
increase the cross section by a factor of 2 to 4 depending on m#± and tan (3. We performed 
a MSSM parameter scan for regions of large cross section assuming a 1 TeV collider and 
a charged Higgs boson too heavy to be pair-produced at such a machine. We find scenar- 
ios with a cross section above 0.1 fb for m#± up to about 750 GeV in the whole scanned 
tan /?-range (1 . . . 50). These scenarios require a sfermion mass scale between 200 GeV and 
600 GeV and a certain amount of mixing in the stop and sbottom sector. The FORTRAN 
program 'eeWH' for the calculation of the MSSM and THDM cross sections including the 
option to perform parameter scans can be obtained from one of the authors 0. 
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Figure 5: Cross section a(e + e~ — > H ± W T ) versus m H ± for different values of tan (3. MSSM 
predictions of the large-stop-mixing and no-stop-mixing scenarios for the unpolarized cross 
section are displayed (thick lines) along with the corresponding sTHDM scenario (thin 
lines) for a collider energy of yfs = 500 GeV and yfs = 1000 GeV. The cross section for 
optimal polarization, P{e~) = —1 and P(e + ) = +1, is shown for the case tan/5 = 15 (thick 
lines with circles). 
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Figure 6: Cross section a(e + e~ — > H ± W T ) versus tan f3 for different values of ra#±. MSSM 
predictions of the large-stop-mixing and no-stop-mixing scenarios for the unpolarized cross 
section are displayed (thick lines) along with the corresponding sTHDM scenario (thin 
lines) for a collider energy of \fs = 500 GeV and ^Js = 1000 GeV. The cross section for 
optimal polarization, P(e~) = —1 and P(e + ) = +1, is shown in each panel for the case 
^_ff± = (thick lines with circles). 
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Figure 7: Cross section for the process e + e — > H ± W T as a function of a common 
superpartner mass scale M for a (a) 500 GeV and (b) 1 TeV e + e~ collider. 
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X t /M f ~ X b /M f ^[TeV] 




M f -[TeV] Xt/Mf X b /M f - 



Figure 8: Two dimensional slices through the scanned parameter space of Scan 1 for 
tan/? = 3, 15 and 30. The colours refer to bins of cross section values: orange: lOfb > 
a > lib, light blue: lfb > a > 10"%, : 10"% > a > 10" 2 fb, cyan: 10" 2 fb > 

a > 10" 3 fb, red: 10" 3 fb > a > 10" 4 fb, green: 10" 4 fb > a > l(r 5 fb, blue: a < 10" 5 fb. 
Points with larger values of a are plotted on top of points with smaller values. 
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Figure 9: Cross section results of Scan 2 versus (a) m H ±, (b) tan/5 and (c) Mi. 
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Figure 10: Two dimensional slices through the scanned parameter space of Scan 2. The 
middle panel (b) shows the same slice as panel (a) but for a variant of Scan 2 which uses 
a simple implementation of the LEP constraint (see text). The quantity X tb in panel (c) 
is defined in Eq. (1231) . The colours refer to bins of cross section values in the same way as 
in Fig. [HI Points with larger values of a are plotted on top of points with smaller values. 
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